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We study the role of the dipolar interaction, correctly accounting for the Dipolar-Induced Reso- 
nance (DIR), in a quasi-one-dimensional system of ultracold bosons. We first show how the DIR 
affects the lowest-energy states of two particles in a harmonic trap. Then, we consider a deep 
optical lattice loaded with ultracold dipolar bosons. We describe this many-body system using 
an atom-dimer extended Bose-Hubbard model. We analyze the impact of the DIR on the phase 
diagram at T = by exact diagonalization of a small-sized system. In particular, the resonance 
strongly modifies the range of parameters for which a mass density wave should occur. 



Introduction. The recent experimental developments 
in the field of ultracold dipolar gases have opened up 
fascinating prospects for the study of systems exhibiting 
Dipolc-Dipolc Interaction (DDI) [l|, Bose-Einstein 
Condensates (BECs) of magnetic atoms have been real- 
ized using Chromium Q , Erbium Q , and Dysprosium Q . 
However, the magnetic moments carried by single atoms 
are small 10 ^b, where /l*b is the Bohr magneton), and 
therefore the effects of the DDI observed with these sys- 
tems have remained perturbative up to now . The re- 
cent realization of the ultracold heteronuclear molecules 
RbK Q and NaK Q, which both carry electric dipole 
moments of the order of one Debye, offer a promising 
route towards stronger DDI effects, but quantum degen- 
eracy still remains to be achieved with these systems. 
Rydberg atoms boast much larger dipole moments and 
will provide an alternative way to explore the DDI, but 
yield challenging experimental problems associated with 
time and length scales @ . 

The DDI is both anisotropic and long-ranged, and 
dipolar gases thus allow for the quantum simulation of 
more general Hamiltonians than the ones accessible with 
non-dipolar neutral particles, whose interaction is de- 
scribed by the standard s-wave interaction [lfj. Fur- 
thermore, the trapping techniques available for ultracold 
atoms makes it possible to control the effective dimen- 
sionality. The use of a lower dimensionality stabilizes 
the system with respect to two-body and many-body 
instabilities caused by the attractive part of the three- 
dimensional DDI. This has prompted detailed studies of 
dipolar systems in 2D and quasi-2D 11 - HJ , bilayer [Hi - 
Il8| . and quasi- id EI-ElI geometries. 

Experiments involving dipolar bosons in optical lat- 
tices have recently been performed both with atomic 
BECs 22 1 and non-condensed dipolar molecules [23j j . 
The simplest theoretical description of such dipolar 
systems in a lattice is provided by an Extended 
Bose-Hubbard Model (EBHM) accounting for nearest- 
neighbour interactions [D]. The ID EBHM has al- 
ready been explored theoretically for generic values of 
the model parameters, revealing the occurrence, beyond 



the standard Mott-Insulator (MI) and superfluid (SF) 
phases, of a Mass Density Wave (MDW) phase [5J |2J| 
and a Haldane Insulator phase [2|| H(§ . 

In this Letter, we study ultracold bosonic dipoles in 
an optical lattice, in the quasi-lD tight-binding regime 
271 . |28| . In order to obtain accurate predictions corre- 
sponding to our specific system, our model should in- 
clude the effect of the Dipolar Induced Resonance (DIR) 
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301 ] . The DIR is a low-energy scattering resonance 
which occurs when the dipole strength is varied. Ac- 
countin g fo r it requires going beyond the single-band 
EBHM [3H. We develop an atom-dimer EBHM, which is 
the simplest possible lattice model reproducing the DIR. 

The scattering and bound-state properties of the DDI 
have been studied numerically for free-space models [32[ 
as well as systems trapped in 3D and 2D optical lattices 
[33I ] . In our quasi-lD geometry, we describe the inter- 
action using an effective potential obtained by averaging 
the transverse degrees of freedom into the corresponding 
harmonic oscillator ground state [ijl, [2^ : 
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where r* = mD 2 /ti 2 is the dipolar length, with D be- 
ing the dipolar strength. The range of this potential is 
determined by the oscillator length l± = (h~/muj±) 1 / 2 
in the strongly-confined directions y and z. The term 
9id = 2h 2 / vnl\ is the strength of the standard s— 
wave contact interaction for a 3D scattering length CI3D 
34| . which can be manipulated using a Fcshbach reso- 
nance [35|]. It competes with the DDI to determine the 
stability and the phase of the system (3|| . For simplicity, 
we assume girj = unless otherwise specified. Note that, 
under this assumption, the potential defined by Eq. (TT]) 
still contains a contact term proportional to r* . 

Two-body physics. The basic building block of our 
many-body lattice Hamiltonian (see Eq. (j4]) below) is 
provided by the solution of the two-body problem in a 
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FIG. 1. (Color online). Left: matrix elements of the two- 
state Hamiltonian instate describing two bosonic dipoles in a 
harmonic trap, as a function of the dipolar length r*. Right: 
corresponding ground-state (red) and excited-state (dashed 
gray) energies. 

single lattice site. Hence, we solve for the ground state of 
two dipolar bosons in a single quasi-lD harmonic well, 
with the axial trapping frequency loq and the correspond- 
ing oscillator length l = (fi/muo) 1 ' 2 . The quasi-lD 
hypothesis requires the ratio 1±/Iq to be small. In this 
two-body problem, the center-of-mass and relative mo- 
tions are decoupled. The relative motion is governed by 
the following Hamiltonian: 

H 2B = p 2 /2m r + m r u$x 2 /2 + Vm(x) , (2) 

where x is the interparticle distance, p is its conjugate 
momentum, and m r = m/2 is the reduced mass. 

Unlike in the case of the contact interaction solved in 
[37J , the Hamiltonian H^b cannot be diagonalized analyt- 
ically. Therefore, we seek its ground state numerically, 
by considering the restriction of H2B onto a subspace 
spanned by a finite number of suitably chosen basis states 
{| (fin)}- Depending on the value of r* , the effective in- 
teraction Vxd supports either no bound state or a single 
one. The bound state is present for large enough values 
of r* , and its entrance coincides with the occurrence of 
the DIR. The 'bare' bound state supported by the at- 
tractive contact part of Vyd(x) plays a key role. The 
corresponding wavefunction is ipBs(x) = \JH exp(— 
where k = r*/(3^), and its cusp at x = cannot be re- 
produced by projecting | ipBS ) onto a subspace spanned 
by any finite number of harmonic oscillator eigenstates 
which are all smooth at x = 0. This implies that the DIR 
physics can only be captured if a wavefunction which has 
a cusp at x = is included in the basis { | <j) n ) } . The 
smallest such basis contains two elements, and we choose 
it to be {| 0o }j I <f>-i )}j where \4>o) is the ground state 
of the ID harmonic oscillator with frequency luq, and 
I 0_i ) oc I ^bs ) — (<po\ ipBS )\4 > o) is a linear combination 
of I tpBS ) and I 4>o ) chosen such that the basis is orthonor- 
mal. Hence, for a given value of the dipolar length r*, 
we replace H2B by the following two-state Hamiltonian: 

„ ( (<M# 2B |<M MH2B\^-l)\ ^ 

il2state - \{4>-i\H 2 B\fo) (<^i|iJ 2B |0-iV • { ) 

The diagonalization of instate (see Fig. [1]) yields approx- 
imations of the ground-state energy E 2 b {t* ) and wave- 
function \^2B(r*))- The energy E 2 b thus obtained is 
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FIG. 2. (Color online). Left: Ground-state energy E2s(r*) 
of the Hamiltonian H2B, for g\n = and l±/lo = 0.8, in- 
cluding 1 (blue), 3 (green), 6 (orange), and 9 (red) harmonic 
oscillator states, without (dashed lines) and with (solid lines) 
the 'bare' bound state | <^>_i ) in the projection basis. Right: 
the four lowest eigenvalues of H2B as a function of r* , for 
QiD = (green) and choosing gin — 2h 2 r* /(3m?5_) to can- 
cel the contact term (red), calculated including four harmonic 
oscillator states and the 'bare' bound state | <f>-i ) in the basis. 

plotted as a function of r* in Fig. fright) for l±/lo = 
0.8. Including a greater number of harmonic oscillator 
states in the basis allows for the calculation of higher- 
energy states (see Fig. [5] (right)). However, as shown on 
Fig. [Deleft), it does not affect the qualitative behavior of 
the ground-state energy E2b(t*) as long as the bound 
state I \ ) is also included in the basis. 

The non-monotonic behaviour of E2b(t*) (Fig. [1] 
(right)) is a signature of the DIR. In particular, the 
ground state energy goes below fruo/2 for r* > r* rit , 
where ^ritAo = 116 for l±/lo = 0.8. Correspondingly, 
the ground-state wavefunction exhibits a cusp, which sig- 
nals a strong contamination from the 'bare' bound state 

I 1 ) . However, the main point for the many-body 

treatment described below is that the population of the 
bound state is very small: the overlap \ (4>-i \ X &2B )| 2 re- 
mains smaller than 0.05 for r* < r* rit . This overlap only 
becomes substantial if the states | 4>o ) and \<f>-i) have 
comparable energies, i.e. for r* /Iq > 6.5 (see Fig.fljleft)). 
The contamination of the wavefunction occurring much 
earlier than the crossing of the bare energy levels can 
be traced back to the shape-resonance nature of the 
DIR [2^], and distinguishes it from Feshbach resonances, 
which usually occur near the crossing of the bare levels 
[3BI ] . More details regarding the scattering properties as- 
sociated with Vid and the DIR will be given elsewhere 

The right-hand plot of Fig. [2] shows the energies of the 
four lowest-energy states of H2B as a function of r* in 
two different situations: (i) the s-wave interaction term 
9id = (green curves) and (ii) the term gm is non-zero 
and chosen such as to cancel the contact term in Eq. ([T]) 
completely [2(J. The r* -dependence of the energy levels 
we predict in these two situations is completely different. 
This will allow for an observation of the DIR using the 
spectroscopic techniques presented, e.g., in [39j . 

Many-body physics. We now consider N dipolar parti- 
cles loaded into a deep quasi-lD optical lattice with unity 
filling factor. We describe this system using a Bosc- 
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Hubbard model [HI Eil extended to include nearest- 
neighbour interactions. We focus on the regime r* < 
r* rit , so that the DIR affects the two-body properties 
even though the number of dimers present in the system 
is extremely small. In order to properly account for the 
resonance, we start from the two-state description intro- 
duced above for the two-body problem (Eq. ([3])). In the 
presence of the coupling between the N different sites of 
the lattice, each of the two states | <j)Q ) an d \(j>-i) yields 
a band and, hence, we introduce an atom-dimer EBHM 
which reads: 

#AD=EJ £ a n i + %rii(ni - 1) - J a (a\a i+1 + he) + Vnin i+1 

(4) 
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+e d m l - Jd{b\b i+ i + he) + Sl{b\aia,i + he)] . 



In Eq. (j4]) , a,i and 6,; are the annihilation operators in the 
site i for atoms and dimers, respectively, and rii = oja, 
and m,i = b\bi are the corresponding number opera- 
tors. The atomic tunnelling coefficient J a is taken from 
4l| . We extract the on-site energy for single atoms and 
dimers, e a and €d, the atomic on-site interaction energy 
U, and the atom-dimer conversion rate f2, from the ma- 
trix elements of i?2stato (Eq. (|3])). The atomic nearest- 
neighbour interaction coefficient V is expressed in terms 
of the effective interaction Vid (Eq. ([T])) and the Gaus- 
sian approximation to the Wannicr function wt (x) local- 
ized on site i. The nearly- vanishing population of the 
dimer band for r* < r cr j t allows for a very crude descrip- 
tion of the dimer dynamics: we neglect atom-dimer and 
dimer-dimer interaction, and we take Jd = J a /10 in our 
numerics. The scaling for Jd follows from the assumption 
that the polarizability of a molecule is twice that of a sin- 
gle atom [4l[. In any case, in the regime we consider, the 
exact value of Jd does not affect our numerical results. 

We focus on the tight-binding regime and we intro- 
duce the harmonic oscillator length Iq characterizing the 
bottom of each lattice well. Like for the two-body prob- 
lem described above, we consider a fixed value of l±/lo, 
which determines the aspect ratio of the system. The 
ground state of the system then depends on two adimcn- 
sional parameters: r* /Iq and V/U. The choice of the 
parameter r* /Iq allows for a direct comparison with the 
two-body physics illustrated in Figs. [1] and [2] In the 
case of the DDI with gm = 0, the ratio V/U does not 
depend on r*. It strongly decays with the lattice depth 
s = Viat/En, where Mat is the intensity of the optical 
lattice and Er is the recoil energy. The validity of the 
harmonic approximation used to determine the single- 
site parameters requires s to be large enough and thus 
imposes an upper bound on V/U. 

For given values of r* /Iq and V/U, we numerically 
calculate the ground state of the EBHM described by 
Eq. (|4]) by exact diagonalization of a 6-atom, 6-well sys- 
tem. In Fig. [H we plot the phase diagram of the system 
as a function of these two variables. The observable is 
the single-particle off-diagonal density matrix element 
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FIG. 3. (Color online). Many-body phase diagrams obtained 
using the single-band (fi = 0, left) and atom-dimer (fi > 0, 
right) EBHMs, performing exact diagonalization on a six- 
atom, six-well system with l±/lo = 0.8. The effective on-site 
interaction U c s < on the right of the vertical dashed line. 



Pi =< a 2 ai >, and it distinguishes the superfluid phase 
(pi 7^ 0) from the insulating phases (p\ = 0). The differ- 
ent insulating phases can subsequently be told apart by 
looking at the projection of the ground-state wavefunc- 
tion onto all accessible configurations. Figure [3] compares 
the phase diagram obtained using the single-band EBHM 
(i.e. taking O = in Eq. (fj|) with the one obtained using 
the atom-dimer EBHM (f2 > 0) for the same range of pa- 
rameters which satisfy the constraints mentioned above. 
Up to a change of units to (U / J, V/J), and in the param- 
eter range for which our model is valid, the single-band 
phase diagram is equivalent to the one presented in [24| , 
and it exhibits two phases: superfluid (SF) and Mott- 
insulator (MI). The two-band phase diagram is quali- 
tatively different. In particular, the MI phase region 
stops at r* = r* rit , and two additional phases appear, 
which are schematically represented in Fig. Upright): a 
(twice degenerate) Mass Density Wave (MDW) phase 
with a spatial period which is twice the lattice period, 
and a 'collapse' phase (which is iV-fold degenerate) where 
all atoms sit in the same well. In our small-sized sys- 
tem, the MI-MDW and MDW-collapse transitions ap- 
pear sharp, in accordance with the first-order character 
expected to accompany a change in discrete symmetry. 
On the contrary, the transitions between the SF phase 
and each insulating phase are smooth, which is compat- 
ible with the predicted Berezinskii-Kosterlitz-Thouless- 
like behaviour expected in ID [24| . 

The phase diagram shown in Figs, fright) andSJlcft) 
can be interpreted in terms of an effective single-band 
EBHM, with the on-site interaction coefficient reproduc- 
ing the two-body ground-state energy: 



H eS ^^ i [e a n i +^n l {ni-l)-J a {a\a l+ i+\Yc)+Vn i n i+ i] , 

(5) 

where U a g(r*) = E2B,{f*) — Tuxiq/I. In the parameter 
range explored on Figs.[3]and|4]deft), the phase diagram 
obtained by exact diagonalization of iJ c ff is very similar 
to the two-band phase diagram. This is due to the de- 
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FIG. 4. (Color online). Left: Zoom-in on the central part 
of the atom-dimer phase diagram (Fig. (right), l±/lo = 
0.8) showing the transitions between the SF, MI, MDW, and 
collapse phases. Right: cartoon representations of the three 
insulating phases for six atoms in six wells. 



tuning A = — U — 2e a between the atom and dimer 
bands being much larger than J7, J a and V 42|. The ef- 
fective single-band model (Eq. ([5])) allows for a direct 
comparison between the phase diagram we present and 
those calculated in terms of the usual EBHM parame- 
ters U/J and V/J. In particular, we find a Haldane-likc 
phase near the upper left corner of our MDW domain, 
which is compatible with the Haldane phase domain re- 
ported in [2(| . A systematic investigation of the Haldane 
phase domain is beyond the scope of the present Let- 
ter and will be carried out on larger-sized systems using 
DMRG [H. 

We also use this effective model to derive quasi- 
analytical approximations of the phase boundaries for 
any number N of particles and sites. For two given con- 
tiguous phases, we calculate the energy of each phase in 
terms of J, U e s, V, and N, using the form of the wave- 
function valid deep within each phase. Equating these 
two energies, we obtain the grey lines shown for N = 6 
on Figs. [3] and |U We now consider more specifically the 
boundary between the SF and the collapse phases, given 
by £ S F-£ colla psc « N{V-2J)-N 2 U cS /2 = 0. The tun- 
nelling term scales linearly with N, whereas the on-site 
interaction term scales with iV 2 . Therefore, for small N, 
the superfiuid phase survives in a region where U a s < 0, 
but the collapse phase is energetically favored for large N. 
This instability favoring the collapse phase corresponds 
to the implosion of a Bose-Einstein condensate with a 
negative scattering length when its size is increased (43| . 

Conclusion. We have studied the effect of the DDI in 
quasi-ID geometries both from the two-body and the 
many-body perspectives. From the two-body point of 
view, we have shown how the DIR affects the depen- 
dence of the energies of the first few excited states of the 
system on the dipolar strength D. Concerning the many- 
body physics, we have pointed out that the DIR makes 
it possible to reach the MDW phase for lower values of 
V/U, and that the system collapses into a single lattice 
site for larger values of D, for which the effective on-site 



interaction U e g becomes negative. The predicted phases 
may be experimentally identified using in-situ imaging 
techniques as well as the recent advances allowing for 
the detection of non-local order 44 1 . 

We studied the effect of the two-body DIR on the 
many-body physics using the minimal EBHM. The DIR 
could also have a strong effect on systems describe d by 
generalized EBHMs such as the one studied in, e.g., [15j . 
It would be interesting to extend this work to 2D geome- 
tries, where the anisotropy of the dipolar interaction is 
expected to play a role. We believe our analysis would 
also be relevant for the understanding of the fermionic 
ID EBHM with repulsive interactions, where the relevant 
phases are the Spin Density Wave (SDW), the Charge 
Density Wave (CDW), and the Bond Order Wave (BOW) 
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